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In this paper we discuss canonical analysis of SO(4, 1) constrained BF theory. The action of this 
theory contains topological terms appended by a term that breaks the gauge symmetry down to the 
Lorentz subgroup of SO(3, 1). The equations of motion of this theory turn out to be the vacuum 
Einstein equations. By solving the B field equations one finds that the action of this theory contains 
not only the standard Einstein-Cartan term, but also the Hoist term proportional to the inverse 
of the Immirzi parameter, as well as a combination of topological invariants. We show that the 
structure of the constraints of a SO(4, 1) constrained BF theory is exactly that of gravity in Hoist 
formulation. We also briefly discuss quantization of the theory. 
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I. INTRODUCTION 

One of the most remarkable developments in general 
relativity of the last decades was Ashtekar's discovery 
that the phase space of gravity can be described with 
the help of a background independent theory of self-dual 
SU(2) connection This discovery became a founda- 
tion of the research program of Loop Quantum Gravity 
0, 0. The original Ashtekar's formulation was gen- 
eralized few years later by Barbero to the case of real 
connections [J], parametrized by a single real number 7, 
called the Immirzi parameter [5j. It turns out that this 
parameter is in fact an additional dimensionless coupling 
constant of the gravitational action, which takes the sym- 
bolic form @ 

S grav = I / e^iW (z*^ 5 + ~ C) > 

(1.1) 

r A = c^P c & a 

e = e t Q /3 7< 5e /J e l/ epe cr . 

In the action above e Q is the tetrad one-form and R a /3 
is the curvature two- form of the Lorentz connection w Q( g, 
where the Lorentz algebra indices a, /3, . . . run from 
to 3. Normally the second term, called the Hoist term, 
regardless of not being a total derivative, does not af- 
fect field equations, because its contribution vanishes on 
shell (for zero torsion) by virtue of the Bianchi identity. 
In spite of this, its presence is not completely innocent: 
it affects canonical structure of the classical theory, and 
quantum theories for different 7 lead to different physi- 
cal predictions (for example the expression for black hole 
entropy calculated in this framework depends on 7 Q). 

It has been noticed in Q that from Wilsonian perspec- 
tive it would be quite unnatural not to append the action 
(|1.1[) with all possible terms that are compatible with 
the field content (e and w) and (local Lorentz and diffeo- 
morphism) symmetries of the theory. It turns out that 



* |rdurka@ift .uni.wroc.pl| 

t jkowalskiglikman@ift.uni.wroc.pl 



there are only three such terms corresponding to three 
topological invariants (Pontryagin, Euler and Nieh-Yan 
classes, see f|4.8|) — f|4. 10[) below). Again, the presence of 
these terms does not influence the classical field equa- 
tion when the constant time slices of the spacetime are 
compact without boundaries. However, they may play 
an important role in quantum theory and/or in the case 
when boundaries are present. In the formulation of Q 
all these terms come with a priori independent coupling 
constants and one wonders if it would be possible to find 
a formulation of the theory so as to organize them in a 
unified way. 

Such a formulation is known for quite some time and 
is dubbed constrained BF theory. The idea that gravity 
can be formulated as a constrained topological BF theory 
has its roots in works of MacDowell and Mansouri [9] and 
of Plebanski [l(| ■ The starting point of the present work 
will be the following action, proposed and discussed in 
[111 (see also [H), 
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is the field strength of the S0(4, 1) (or S0(3, 2)) connec- 
tion A^ IJ , while B^ IJ is a two-form field valued in the 
algebra of the same gauge group. The capital Latin in- 
dices I,J,K,... are the algebra ones and run from to 
4, when the Lorentz subalgebra of the gauge algebra is 
labeled by Greek indices from the beginning of the alpha- 
bet a, /3, 7, . . . running from to 3. We will decompose 
them into timelike and spacelike a,b,c, . . .. Below, in 
the course of Hamiltonian analysis we also decompose 
the spacetime indices /i, v into time and space denoting 
the space indices by letters from the middle of the Latin 
alphabet i,j, k, . . .. 

As we will show in the next section, the theory defined 
by the action (II. 2|) is equivalent to Einstein-Cartan the- 
ory with action accompanied with the Hoist term and 
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the topological terms described above. The six coupling 
constants of Q are then replaced by two dimensionless 
couplings a and (i of (|1.2j) and one dimensionful scale I. 
In the Sec. II we will discuss the canonical formulation 
of this theory, while in Sees. IV and V we will show 
how these constraints can be simplified and recast into 
the form proposed by Hoist. In the final section we will 
make some comments concerning perturbative quantiza- 
tion of the theory around Kodama state. 



II. GRAVITY AS A CONSTRAINED BF 
THEORY 

In this section we will recall some properties of the 
action (jl.2j) . It has been shown in [TTJ] that this action 
is equivalent to the standard action of Einstcin-Cartan 
gravity. To see this one first decomposes the connection 
A p IJ into tetrad and Lorentz connection 



(1:4 



A a P 



a/3 



(2.1) 



with £ being a length scale, necessary for dimensional 
reasons since the connection on the left hand side has 
the dimension of inverse length, while tetrad is dimen- 
sionles^j] associated with the cosmological constant 



1 

£2 



A 
~3 



Then one solves equations of motion for B and substi- 
tutes the result back into the action. As a result one 
finds Einstein action appended with a number of topo- 
logical invariants. To find its canonical form one has to 
associate the dimensionless coupling constants a and /3 
of (|1.2[) with the physical ones: Newton's constant G, the 
cosmological constant A, and the Immirzi parameter 7: 



GA 



a = 



3 (l + T 2 )' 



D-^jrl^, 7-5. (2.2) 
3 (I+7 2 ) a 



In (|2.3p T>^ is the covariant derivative defined by connec- 
tion A, so that 

(D p B vp ) IJ = d p B vp IJ + A^k B up KJ + A^k B vp IK . 

The theory defined by (11.21) for non-zero a breaks the 
original de Sitter SO(4, 1) gauge symmetry down to 
Lorentz SO(3, 1). It is, therefore, convenient to decom- 
pose the covariant derivative T> p into Lorentz so(3, 1) 
and translational parts, and to use the Lorentz covariant 
derivative defined by Lorentz connection u) (|2.1I) , to wit 



(2.5) 

(P* B vp )^ = (V» B vp ) aA - 1 e pfs B vp a P , (2.6) 
where 



{ViB^Y? = d»B v 



(2.7) 

with an obvious generalization for another Lorentz ten- 
sors. Using this decomposition we rewrite the field equa- 
tions d%U) , as 

ciivp „ ( jy* Bv * _ 1 e a BvpPi + 1 e ^ Bu rt\ = } 

(2.8) 

< (m B up ai - j e p0 B vp ^\ = , (2.9) 



V - (3 B^ - - e a ^ s B pv jS = , (2.10) 



f3B^ a4 = 0. 



(2-11) 



Notice that the curvature in f|2 . 10[) is the sum of Riemann 
tensor of u and the cosmological curvature 



Instead of repeating this derivation here, let us show 
that field equations resulting from the action (ll.2|) are the 
standard vacuum Einstein equations. The field equations 
read 
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(2.4) 
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while that in (|2.1ip is just the torsion 



Solving (|2~TQ1 and ([2~TT| for B we find 



(2.12) 



V 4 = \<Pt ^ a ~ v$ O = J V ■ ( 2 - 13 ) 



where 



1 In our approach all generators of the gauge algebra are dimen- 
sionless. Alternatively, one can use dimensionful generators of 
the translational part of the algebra (as it is usually done when 
one wants eventually to make the algebra contraction). Then 
momentum generators have canonical dimension of inverse length 
and I shows up in the algebra as well. 



with 8"$ = 5%8$ - 8^6 f . The tensor M is a sum of 
Lorentz invariant tensors and, therefore, its covariant 
derivative T>f, vanishes. 
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Substituting (|2T4| into (JUJ and using Bianchi iden- 
tity for Riemann curvature one can check that the result- 
ing equation forces torsion T^ u a = iF^" 4 to vanis 
Using this it is easy to see that (|2.9p is equivalent to 
Einstein equations with cosmological constant A = 3/£ 2 . 
This completes the proof that field equations following 
from the action (|1.2jl reproduce the standard Einstein 
equations. 

It should be noticed that when the coupling constant 
a = the theory becomes topological, so that the last 
term in the action (|1.2j) that explicitly breaks the gauge 
symmetry from the topological SO(4, 1) down to physi- 
cal S0(3, 1) carries all the information about dynamical 
local degrees of freedom of gravity. As we will see below 
this fact is clearly reflected in the structure of constraints 
algebra. 

III. CANONICAL ANALYSIS 

In the first step of canonical analysis of the constrained 
BF theory defined by (|1.2[> let us decompose the curva- 
ture F lll , IJ into electric and magnetic parts 

F^," ^(F Qi IJ ,F iJ IJ ) (3.1) 

with 

F 0l IJ = A IJ - d.Ao" + Vx A kj - Vk A kj 

= A/ J - 2W J (3.2) 

where the dot denotes the time derivative, T>i is the co- 
variant derivative for the connection Ai IJ , and 

F tj IJ = ft V 7 + A iK Aj KJ -mj. (3.3) 

As usual the zero component of the connection becomes 
a Lagrange multiplier for Gauss law. Further we decom- 
pose B field into 

B^ IJ -> (B 0l IJ = B/ J , V lIJ = 2e ijk B jk IJ ) . (3.4) 

As we will see shortly, V lIJ turn out to be momenta as- 
sociated with spacial components of gauge field A, while 
the remaining components of B play a role of Lagrange 
multipliers. 

Using these definitions and integrating by parts we can 
rewrite the action as follows 

5 = I dtL, (3.5) 



L = / (fx^uAi" + B/^u + A IJ n LJ ) . (3.6) 



It is clear that and Ao IJ are Lagrange multipliers 

enforcing the constraints lT/j and II/j, which explicitly 
read: 

Hu(x) = {v i v i ) I1 (x) 

= (diV^j + A u K V % K j + A U K V % IK ) (x) (3.7) 

which is the Gauss law for S0(4, 1) invariance (see be- 
low), and 

U'uix) = (2e^ k F jkIJ - P Vh-« 2 e IJKLA V iKL ) (x) 

(3.8) 

The Poisson bracket of the theory is 

WW^klM} = \&ix- V) 4 &KL • (3-9) 

(The factor 1/2 results from the fact that the canonical 
momentum associated with A defined as 5L/SA is 2V, 
not P.) The Lagrangian (13.6[) contains just the standard 
(pq) kinetic term appended with a combination of con- 
straints, reflecting the manifestation of diffeomorphism 
invariance of the action (|1.2I) that we have started with. 
It is worth noticing that prior to taking care of the con- 
straints the dimension of phase space of the system is 
2x3x10 = 60 at each space point. As we will see the 
dimension of the physical phase space is going to be 4, 
as it should be. 

The Poisson brackets of the constraints can be straight- 
forwardly computed and read 

{U ij (x),IVkl{v)} = 5{x - y) (n IL ILj K (x) 

- Vjl^ik{x) - vikUjUx) + ru K R IL (xj) « (3.10) 

which means that 11/ j form a representation of the gauge 
group 50(4, 1) of the unconstrained theory (a = 0), as 
expected. Further 

{Wu{x\W KL {y)} = 2ae^ k 5{x - y) (e KLIPi A k p j(x) 

- £KL.jPiA k F ' i{x) + eijKPiA k p L {x) - ei.jLPiA k p k{x)\ 

(3.11) 

and 

PwW, nVi(y)} = ~6{x-y) (e KLIPi V lP j(x) 
-£KLjPiP lP + ^S(x - y) (jlitejKMNA 
—"HJLtlKMNA — VlK^JLMNi + VJK^ILMNij'P lMN (x) 

+ o ™ V) {vil^A 1 jk{x) - rjj L II l iK(x) 



2 To prove this one has to assume invertibility of the tetrad. 



- Vik^Mx) + VJKtfidx)) . (3.12) 
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It is worth noticing that in the topological limit a — all 
the constraints are first class. This observation leads to 
the following, apparent puzzle. Namely, as we said above 
the kinematical phase space is 60 dimensional. On the 
other hand for a — we have 10+30 first class constraints 
that remove from this phase space 80 degrees of freedom. 
How is this possible? To answer this let us notice that not 
all the constraints are independent. Indeed taking the 
covariant divergence of the H l u constraint and making 
use of the Bianchi identity we see that 

{ViU^j = -pn u (3.13) 

and thus the set of constraints is reducible. It follows 
that we have only 30 independent first class constraints 
IP/j, which remove exactly 60 dimensions from the phase 
space, as it should be since the theory with a = is 
topological. 

Returning to the case a^0 we notice that the action 
(|1.2[) is invariant under local gauge transformations that 
belong to the Lorentz subgroup SO(3, 1) of the initial de 
Sitter group SO(4, 1)0. It follows that it is natural to 
expect that one can simplify the algebra of constraints 
(|3.10[) - (I3.12[) if one decomposes the constraints into that 
belonging to the Lorentz and the translational parts of 
the algebra. From (I3.7P we get 

n a4 (x) = U a (x) = {VtV i ) aA (x) - je/(x)V\ (x) « 

(3.14) 

U t ,p(x) = (VfP i ) ap (x) 

- je la (x)V l pi (x) + je lP {x)V l ai {x) « (3.15) 
while from ([3T5]l 

ir Q4 (x) ee <S>\{x) = (2e^' fc F jkcA {x) - PVUi*)) « 

(3.16) 

W a p{x) EE $ l Q/3 (x) = 

= (2e« fe F ikaP {x) -Wipix) - \ e aM5 V^ s {x)) « 

(3.17) 

One then finds that the algebra of constraints (|3.11l) . 
p,12j) simplifies a lot, and the only brackets that do not 
vanish weakly are 

{Il a (a!), $V(2/)} « -f S(x - y) e lSap V lp i(x) (3.18) 
{U a (x), <S>\(y)} « -~ S(x - y) e aipa V^(x) (3.19) 



In what follows we restrict ourself to the positive cosmological 
constant case; the negative cosmological constant and the Anti 
de Sitter group SO(3, 2) can be analyzed analogously. 



{$ l Q (a;),$V(y)} ~2ae ljk S(x -y)e al 3 7S A k s i (x) = 

^e^ k 6(x-y)e aMS e k 5 (x) (3.20) 

Now we can turn to the next step of canonical analysis, 
i.e., to checking if there are any tertiary constraints. The 
Hamiltonian, being a combination of constraints reads 

H = -2A a n a -A^n aP -2B l a ^ a -B t a ^\p (3.21) 

It follows from (|3.18fl3.20|l that we have to satisfy the 
following conditions to ensure that the constraints are 
preserved by time evolution, generated by hamiltonian 

Ha = f (BS V^ 5 + Bi*> P l \) e aPlS w (3.22) 
= ~a (-^ k Bp e k s - ^V^ 5 )e aPlS « (3.23) 

= -oi (je l]k Bf e k s + A^V u A )e aM& « (3.24) 

These equations can be solved for Lagrange multipliers 
(we have 34 equations for 34 unknowns , B^ 1 , A 1 
with arbitrary coefficients) and thus there are no tertiary 
constraints. 

Notice however that there is an ambiguity in Dirac 
procedure in the case of diff-invariant systems, i.e., such 
that hamiltonian is a combination of constraints. The 
usual approach is to check if one can solve the vanishing 
of time derivative of the constraints condition for La- 
grange multipliers, as we did above. But this is, clearly, 
not a general solution of these conditions. In general one 
may look for the solutions with arbitrary values of the 
Lagrange multipliers, but instead restricting the phase 
space (for example if we impose the condition that all 
the Lagrange multipliers in (|3.22p ~ (|3.24p are arbitrary 
there would be additional constraints saying that com- 
ponents of tetrad and momenta are to be equal zero.) 
Notice that this problem does not arise in the case of 
the hamiltonian not being weakly zero, because then the 
resulting equations pertaining to the time invariance of 
the constraints are non- homogeneous. Thus the proce- 
dure that is usually employed does not seem to provide a 
complete characterization of the phase space, but we will 
adopt it here, leaving the discussion of this subtle point 
to the future work. 



IV. SIMPLIFYING THE CONSTRAINTS 

The aim of this section is to rewrite the system of con- 
straints (|3.10p - (|3.12l) in a form that makes comparison 
with constraints of General Relativity with Hoist term, 
discussed in [6j. In what follows we will borrow some 
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ideas from the paper of Perez and Rezende Q. (Similar 
ideas, albeit in more restricted setting, were discussed, 
e.g., in [H and 0].) 

In the first step let us rearrange the constraints (|3.10p - 
(13 . 12[) to write them in the following form 



=P' i . - M- n1 S P.. . c ^ 



-a/3 - I'afi - M a/3 Fjk 7< 5 £ 

2 



n Q/3 = -^e^V^K^e^eks) » 



(4.1) 
(4.2) 
(4.3) 



2q 



(a 2 + /3 2 )£ 3 



e tf *ea^e/e/e fc a »0 (4.4) 



Recall that the coupling constant a and f3 satisfy the 
identity a /(a 2 + /3 2 ) — ( 2 /G, while the operators M and 
if are defined to be 



Ma ^ S= (a 2 +/? 2 ) V7 ^ 



(7^-^), (4.5) 



7<5 



(a 2 +^ 2 ) V 7 ^ 



+e«V)- (4-6) 



Also recall that the action (|1.2I) . after solving for i? and 
expressing the resulting action in terms of the SO(3, 1)- 
connection ui and tetrad e, has the form [ll[ 

S — — f f al3 '* s (R Re p x - — e e ae e x\e iiVpa 



7 G 



2GA 4GA 



(4.7) 



One immediately recognizes here the standard gravita- 
tional action in the first line, and the Hoist term, whose 
strength is governed by the Immirzi parameter 7 = /3/a 
in the second. The last three terms are proportional to 
topological invariants (Nieh-Yan, Ponryagin, and Euler): 



NY A = J (T M „ a T° - 2 a pe?eP) t^ pa , (4.8) 
Pa = J R^apB^eT^, (4.9) 
E± = J R^ a pR P ^se a ^ S e^ . (4.10) 

As we will show, in the case when the constant time 
surface is without boundaries <9£ = 0, the topological 
terms play the role of the generating functional for canon- 
ical transformations, which simplify the constraints con- 
siderably Q . The key observation is that Pontryagin and 
Nieh-Yan invariants can be expressed as total derivatives 

NY A = 4 Jd^evcPZe?) (4.11) 
P 4 = 4 fd p (u v at d p ^ + |w„ ab u; e u?) e^ pff (4.12) 



The same holds for Euler class. However in this case one 
has to make use of self and anti-self dual combinations of 
Lorentz connection 

1 



(4.13) 



and curvature (see e.g., [15[) 
1 



- - ( R a P m - f a V W s \ 



(4.14) 



It can be checked that both Pontryagin and Euler class 
can be rewritten with the help of ± -RJ^ as follows 

Pi = J e^ p (+R^+R pa afj + -Rf v -R pa afi ) (4.15) 



E 4 = 2i J e^P(+R a £+R paaP --Rf v -R paaP ) (4.16) 
Introducing 

C» = (w„ d^f + lu v ab cj p a c cjf) (4.17) 

we write Pontryagin and Euler classes as total derivatives 

P 4 = 4 J (d p C p (+oj) + d p C p {-u)) (4.18) 

E 4 = 8iJ (d p C p {+uj) - a M C(-w)) (4.19) 

Therefore the topological part of action (|4.7p takes the 
form 



It is worth noticing that in spite of the presence of the 
imaginary i here, the action St is real (for real 7.) 

For constant time surfaces, being a manifold without 
boundary (<9£ = 0), all total spacial derivatives terms 
drop out and only the ones with total time derivative 
survive 



S T = J d W{e,uj), 

where W{uj, e) is a functional of torsion and self and anti- 
self dual Chern-Simons forms Ccs = C a 



(4.21) 



2a 



(a 2 + /3 2 ) 



((7 - i)C C s{ + u) + (7 + ^cs(^)) . 
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Having the functional W we can make canonical trans- 
formation, which defines new momenta J 21 *, £? % ah of the 
tetrad e and the connection oj, respectively 

K = KHK> W(u, e)}, ^ = V^HKp, W(u, e)} 

(4.22) 

with 

{e?,&i} = \l8i5<$ and {u>? = ±% . 

(4.23) 

Since the variations of the functional W{uj, e) are 

1 x ,4 

2 ^ = % 7^ ^ --ppe ja e kp (4.24) 



2 8e? £B Jk 



(4.25) 



we find that the resulting constraints, expressed in terms 
of new momenta (|4. 22[) take the form 



K =^«0, (4.26) 
^ = ^ - i V e, 7 e k s e^ k « (4.27) 



IW = ^ e« fe # a/5 ^ (e J7 e^) w (4.28) 

(4.29) 



This form of constraints will be our starting point in 
checking equivalence with the ones proposed by Hoist 
Q, which will prove, in turn, that they describe General 
Relativity, as expected. To establish this equivalence we 
will have to fix the time gauge. We will turn to this 
problem in the next section. 

Before closing this section let us make an important re- 
mark. The considerable simplification of the constraints 
relies heavily on the fact that the constant time sur- 
faces are manifolds without boundary. In the case when 
boundaries are present the analysis of the constraints be- 
comes much more involved. We will address this issue in 
the forthcoming paper. 



V. TIME GAUGE 

In order to make contact with the Hamiltonian analysis 
of Hoist, we have to fix the gauge so as to remove the 
time component of the tetrad and then to relate momenta 
associated with Lorentz connection with an appropriate 
combination of the remaining tetrad components. 

To this end, let us introduce the gauge condition which 
must be added to the list of constraints 



(|4~26|) . (|4~27l) take the form 
2a 



(5.2) 



- pjf+p^ ft ^ o e kd « (5.4) 
where we have used the convention e 0abc — e abc and 

^Oabc ^abc- 

Combining the last two equations we find constraints 
for generalized self and anti self-dual parts of & 



-Pi «o 



(5.5) 



4a 



where we define 



(5.7) 



It can be easily checked that ± P are momenta associated 
with generalized (anti) self-dual combinations of Lorentz 
connection (which for 7 — ±z become usual self and anti 
self dual ones) 

±„.,a Oa _i_ \ ^"-bc (r Q\ 

w i — w i ± 7T e UJ lbc (O.OJ 
2 7 

with the Poisson brackets being 

±P* } = , ±P^} = 5j6 a d . (5.9) 



Let us now turn to the constraint H a p (|4.28l) . Decom- 
posing it into components we find 

U "b = ( a 2 +^2)£2 ^ ( - 2£abc d e d e k + 



2 \ 
H — {di{e ja e k b) +w iQ c e jc e kb - ui lb c e 3 c e ka )J (5.10) 



n 0Q 



-2a 



ijk 



{a 2 + ii 2 )t 2 

(e abc {d i {e)e%)+2u\ d e jd el) + -w 4 Q b e, b e ka ) (5.11) 



7 , 

Taking the combination Ho a ± ^abc^- we get 



4a 



(a 2 + /3 2 )£ 2 



ijk 



and 







(5.1) 



In this gauge the constraints w can be removed by 
turning to the Dirac bracket so the remaining constraints 



4a 



{a 2 + P 2 )P 



ijk 



, 1-7% b 

{ ) u i0 e ja e kb 



e abc {d i (e b e%)+2cj bd e jd e%))^0. (5.13) 
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Prom these two equations it follows thal[j| 

^ 06 e JQ e fc fe 6^«0 (5.14) 
(d^e) e%) + 2u\ d e d e c fe )e^ e abc a (5.15) 

which expressed in terms of the new variables 
1 



^ a = 5 ( + < + -<) , = 



) 

(5.16) 

take the form of the Gauss and the boost constraints 

G a = { + w\ + -wt) e 2 a e k b e l ' k « , (5.17) 



B a = e c k ) e a bc - 1 { + w\ ~ w\) e 3 a e k b)> 



ijk 



We can handle the scalar part of (|4.29[) 
S 



(5.18) 



(-elR jk0c -e abc e^Fj k c )e i ^^O 



(a 2 + P 2 )t V 7 
similarly, obtaining as a result the expression 



(5.19) 



\ 7 



)e dbc d 3 { + Wkd ) + (L_f.y^ dj[ - Wkd) 



>+7 2 W.,+„.. c <■■>> A ,, ,, 



+ (l+ 7 2 )+ Wj & -<-(i±^-)+^+<-(^)- 



a ijfc 



(a 2 + /3 2 )^ 



0. 



As for the vector part of (|4.29l) 



V a = 



we find 



e\ R jk ab + 2e ab c e\ R jk 0c ) e» k « , 

(5.20) 



1 + 7 \,+ - - + x 

2 ^ J { Wj b W kc + Wjb^Wkc) 



c f ijk ^ ~ Q 



e,- e 



G 



It should be noted that in the case 7 2 = — 1 all terms 
containing + m; a cancel. Notice also that with the help of 
the Gauss constraint the vector constraint can be reduced 
to the form 



2e ai)C 3- w k 



i-tV 



w jb w kc 



e? e ijk « . 



The form of the constraints that we have obtained so 
far is still not the final one. At some point we will have to 
get rid of the constraint + F' a and eliminate the de- 
pendence on + Wj b of all the remaining constraints. How- 
ever in order to be able to do that we must simplify the 
form of the boost constraint (|5. 18|) . To see how this can 
be done we multiply the constraint (|5.5j) by tetrad and 
decompose the resulting constraint C a b into symmetric 
and antisymmetric parts 

C (q6) = + F, (a ei } , C [ab] = +P, [tt 4] . (5.21) 

Let us now calculate the Poisson bracket of C ab with the 
scalar constraint 

{C a b, S} = 



i b{^d k e ia -Y + w c k --wt)e d i e a Jj . (5.22) 



(i + 7 2 K c » fe( 

The bracket of the antisymmetric part C\ a b] gives exactly 
the boost constrain^ (|5 . 18[) . However, the bracket of the 
symmetric part C( q m leads to the secondary constraint 

B ab = 2e ljk (e j bd k e ia + e ja d k e lb ) 

-j( + w c k - ~w c k )e d {e jb e acd + e ja e bcd ) rjO. (5.23) 

Clearly, this constraint would arise if we impose the re- 
quirement that all the constraints are to be preserved in 
the time evolution. Therefore one has to add B ab to the 
set of constraints of the theory. But then the road sud- 
denly becomes sunny. It suffices to note that the boost 
constraints B a and the newly derived constraints B ab are 
just the antisymmetric and symmetric parts of the simple 
constraint 



2r? 



7( + < 



0. 



(5.24) 



where Tf is a (unique) solution of the Cartan first struc- 
tural equation 

(d {i e j]b ) + eb cd Y c {i e d ^ =0. (5.25) 

Using (pT2"4")l and (jSTSj) we get rid of both +P' a and 
+ uij bl replacing + Wj b in all the remaining constraints 
with the solution of (|5.24[) . Similarly using (|5.2j) and 
(|5.6p we can identify the momentum of ~Wj b with 

4a Jjk , u p b p c 

(a 2 +p 2 )l 2 e *bce 3 e k . 

What remains are therefore 3 Gauss, 3 vector and 1 scalar 
constraints, all of them first class, constraining the 18- 
dimensional phase space of ~Wj b and its momenta. Thus 
the dimension of physical phase space is 18 — 14 = 4 as 
it should. Of course, the final set of constraints we have 
obtained has exactly the form of the constraints describ- 
ing gravity, cf. @. This completes our analysis of the 
canonical structure of SO(4, 1) constrained BF theory. 



4 This is obvious for 7 2 ^ -1. For 7 2 = -1 eq. 15.12) is identically 5 Instead of C[ab] = + P; [a e* just take the expression C ab e abc = 
satisfied, but then, since the connection is real, the real and 2 
imaginary parts of J57i3t lead to l5A4t . (Q5) . +P« a 4 so {°ab t abc , S} = ' B c . 
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VI. COMMENTS ON QUANTIZATION 

Let us conclude this paper with some comments con- 
cerning quantization. Clearly, one can take the first class 
Gauss, vector, and scalar constraints as a starting point 
in construction of the quantum theory, as it is done in 
Loop Quantum Gravity |2|], [3|. However, the structure 
of constraints of the original theory opens another pos- 
sibility of devising a perturbative expansion in parame- 
ter a around topological vacuum. Here we will describe 
briefly this perturbative theory leaving details to a sepa- 
rate publication. 

Our starting point will be the set of constraints (|4.ip ~ 
(|4.4[) . Consider now the canonical transformation (|4.22l) . 
Its quantum counterpart can be easily found. To see 
how, take the mechanical model in which one makes the 
transformation (see PH) 

Pi ->Pi =Pi + {Pi,f(0)} 

so that quantum mechanically we have 

Pi +i\pi,f(q)] ■ 

If we represent pi = id/dq l thenp- = id/ dq 1 — df{q) / dq l . 
Therefore if we decompose the wave function ip(q) — 
exp(-if(q))ip'(q) then 

p'ip{q) = exp(-if(q))p'ip'(q) , 

which means that we just have to multiply the wave 
function with the phase exp(— if(q)) and then use the 
standard representation of the new momenta p' as the 
derivatives over positions. In the case at hands (|4.22|) . 
it is therefore sufficient to multiply the wave function by 
the prefactor exp (— iW(e, cj)) where W(e,u) is given by 
(|4.22p , and replace all the momenta V with the new ones 
Then we can just use the constraints ()4.26|) - (|4.29|) . 
When a = these constraints reduce to the first class 
set 



0, ^ 



0. 



(6.1) 



The wave function annihilated by them is just a con- 
stant, and thus the full physical wave function is a phase 
exp (— iW(e,u>)). Clearly, and not surprisingly, in this 
case the wave function is the Kodama state [17j (strictly 
speaking this is the Kodama state for SO(4, 1) multiplied 
by the phase proportional to Euler class of a constant 
time manifold.) Notice that here this state is delta func- 
tion normalizable, because all our constraints are real (cf. 
[l8j). The simplicity of the zeroth order (in a) solution 
reflects the fact that to this order the theory is topologi- 
cal. 

Let us now turn to devising the a perturbative the- 
ory. The constraints (|4.26[) - (|4.29[) all have the form 
$ = $(o) (for the last two $(°) = 0). We also 

expand the wave function in the series in a, to wit 

* = * (0) +a* (1) + ... . (6.2) 



The problem we are facing now is that for non-zero a 
the constraints are no longer first class and therefore we 
need a nonstandard procedure to handle them. One pos- 
sibility would be Gupta-Bleuler quantization (l9j . but 
the required procedure of splitting the constraints into 
holomorphic and anti-holomorphic parts is technically 
complex and, presumably leads to explicit breaking of 
Lorentz covariance (see [20| for discussion in a similar 
context.) Another possibility would be to make use of 
the master constraint program [2l|, [HI, [Hj], and [20| . 
but this is again technically involved. 

Instead we adopt the definition of physical wave func- 
tion ^ such that the matrix elements of all the constraints 



(*|$|*) =0, 



(6.3) 



which is a weakened version of Gupta-Bleuler scheme. 
It should be stressed that the expression (|6.3I) is formal, 
because to make the precise sense of it we must specify 
the inner product in the Hilbert space of states. 

Now we use (|6.3p to define the perturbative theory in 
a. In the zeroth order we have 



^(0) $(o) #(0)\ = o. 



while in the first order in a we find 

$(°) ^(0)^ + ^(0) $ (o) ^(l) 



(6.4) 



+ 



$(i) 



0. 



(6.5) 



Inspecting (|4.26p ~ (|4.29l) we find that it follows from ( 
(I6.5P that the zeroth order wave function has to satisfy 
the following four conditions 



0= (*(°) 
= 



5ef(x) 
S 

"WW 



e* k Kef 6 ef F, 



jk 7(5 



(6.6) 

(6.7) 

(6.8) 
(6.9) 



Knowing one can turn to the remaining first order 
equation, resulting from ([4.271) . along with some of the 
second order ones, to find and then go to the next 

order analysis. We stop the discussion at this point leav- 
ing the details to another paper. 



VII. CONCLUSIONS 

In this paper we have performed the canonical analy- 
sis of the constrained SO (4, 1) BF theory. This analysis, 
although quite involved, seems to be significantly sim- 
pler than the analogous one of Plebanski theory reported 



9 



in [24[, leading however to the slightly more general ef- 
fective description of the dynamical degrees of freedom 
provided by Hoist constraints that include Immirzi pa- 
rameter. This suggests that it might be not only simpler, 
but also more natural to consider spin foam model asso- 
ciated with this particular formulation of gravity. Unfor- 
tunately, not much work has been done till now on the 
SO(4, 1) spin foam models, which would require to han- 
dle somehow not only the quadratic B field term, but 
also the representation theory of SO(4, 1) group, which 
is more complicated than the one of SU(2) group, usually 
used in the spin foam context. 
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